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Motivated from the experimental observation of metal induced gap states (MIGS) at insula- 
tor/metal interfaces by Kiguchi et al. [Phys. Rev. Lett. 90, 196803 (2003)], we have theoretically 
investigated the electronic properties of MIGS at interfaces between various alkali halides and a 
metal represented by a jellium with the first-principles density functional method. We have found 
that, on top of the usual evanescent state, MIGS generally have a long tail on halogen sites with a 
Pz-like character, whose penetration depth (A) is as large as half the lattice constant of bulk alkali 
halides. This implies that A, while little dependent on the carrier density in the jellium, is dominated 
by the lattice constant (hence by energy gap) of the alkali halide, where Alif < Alici < Alii- We also 
propose a possibility of the MIGS working favorably for the exciton-mediated superconductivity. 

PACS numbers: 73.20.-r, 73.40. Ns, 71.15Mb 



I. INTRODUCTION 

Heterointerfaces, namely solid-solid interfaces between 
very dissimilar materials, have been an issue of great 
interest in condensed matter physics, where various 
fascinating possibilities may be considered. Espe- 
cially, the electronic structure of metal/insulator or 
metal/semiconductor interfaces has a long history of in- 
tensive studies dating back to the 1930's^'^. 

One of the most crucial factors which govern the elec- 
tronic structure of heterointerfaces is the metal-induced 
gap states (MIGS), whose eigenenergies lie in the band 
gap of the semiconductor (or insulator). While MIGS 
spread well into the metal side of the interface, they decay 
exponentially into the semiconductor. The decay length 
is an order of a few A for typical semiconductors, where 
MIGS have been believed to determine the Schottky bar- 
rier height'^"^. 

If we have an insulator/metal interface rather than 
a semiconductor/metal interface, on the other hand, it 
has generally been believed that the decay length of the 
MIGS into the insulator is negligible. This view has been 
based on quite a plausible model calculation, where the 
insulator/metal is represented by a simple tight-binding 
model/jellium^. For a typical ionic insulator such as 
LiCl, the band gap (~ 0(10) eV) is an order of mag- 
nitude greater than the transfer energy (^ 0(1) eV) be- 
tween the cations and anions, and the decay length is 
estimated to be smaller than 1 A, which would suggest 
that MIGS should be irrelevant at insulator/metal inter- 
faces. 

Experimentally, it is difficult to detect MIGS at the 
heterointerface between wide gap insulators and metal, 
because clean interfaces are difficult to be fabricated. 
While MIGS have been reported in MgO/Cu(222)io or 



MgO/Ag(001)^^, the question of whether the observed 
gap states are formed by local chemical bonding re- 
mained, since MgO(222) is a polar surface so that a 
strong hybridization between the O 2p band and the Ag 
5sp band is expected for MgO/Ag^'^. 

Recently, Kiguchi et af^ have succeeded in observing 
MIGS at LiCl/Gu and LiGl/Ag interfaces. They have 
grown a LiGl film on a metal substrate in a layer-by-layer 
fashion, and exploiting the near edge x-ray absorption 
fine structure (NEXAFS) which is an element-selective 
method, they have observed an extra peak below the 
bulk edge onset for LiGl film, which was interpreted as 
an evidence for MIGS. The reason why the detection 
eluded earlier studies is that, in conventional methods 
such as ultraviolet photoemission spectroscopy (UPS), in- 
verse photoemission spectroscopy, or electron energy loss 
spectroscopy (EELS), the signal from the interface is ob- 
scured by significant contributions from the substrate. 
Kiguchi et al. have estimated the decay length of MIGS, 
which are 2.6 ± 0.3 A for LiGl/Cu(001) and 2.9 ± 0.7 
A for LiCl/Ag(001), values much greater than those ex- 
pected in simple models^. An ab initio electronic struc- 
ture calculation for LiGl/metal system has supported the 
existence of the MIGS that is strongly localized at the 
interface^^. 

The purpose of the present paper is to study the 
electronic structure of interfaces between various alkali 
halides and a metal by means of first principles calcu- 
lation, and to search for novel possibilities specific to 
metal/insulator interfaces. We shall show that (i) MIGS 
have a tail on halogen sites with a Pz-like character, 
whose penetration depth (A) is as large as half the lat- 
tice constant of the bulk alkali halide, which are totally 
consistent with the experimental result^^, (ii) while A is 
insensitive to the carrier density in the metal, the pene- 
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tration depth is dominated by the energy gap of the alkah 
halide. SpecificaUy, the band gaps are ii^LiF > £^LiCi > 
i?Lii, and the the penetration depth has the expected in- 
equahty, Xuf < Alici < Alii. 

We can consider various possibihties for many-body ef- 
fects unique to heterointerfaces. There is in fact a long 
history for those, where the notable ones include metal- 
insulator transition^'' and superconductivity'^^'^^'^®'^^. In 
the context of the latter, here we propose that large den- 
sity of states arising from the MIGS inside the insulator is 
a good news for the superconductivity due to the exciton 
mechanism. The history of the exciton mechanism of su- 
perconductivity dates back to the proposals by Little^^, 
Ginzburg^^'^*, and Bardeen et af^. For this mechanism, 
the coexistence, in real space, of dispersionless excitons 
and metallic carriers is required^^. While this is difficult 
to achieve in general, if the density of states due to MIGS 
inside the insulator is large enough, we may envisage that 
the insulator/metal interface such as those discussed here 
may provide a possible ground for superconductivity. We 
shall discuss this toward the end of the paper. 
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FIG. 1. Front(a,b) and top(c) views of the atomic config- 
uration considered in the present study. In a(b) one (two) 
layer(s) of alkali halide sandwich a metal region modeled by 
a jellium, where solid (open) circles denote cation (anion) 
atoms while shaded regions the jellium. Unit cells are indi- 
cated. The lattice constant along x, y directions is denoted 
as ax = ay = a/\/2, where a is the lattice constant defined 
for the bulk alkali halide with NaCl structure. The interlayer 
distance is denoted as a^ = a/2. 



II. METHOD 



The electronic structure of MIGS is studied with a 
first-principles calculation based on the density func- 
tional theory. We adopt the exchange-correlation func- 
tional introduced by Perdew, Burke and Wang'^^ and 
ultra-soft pseudo-potentials^"'^^ in a separable form. We 
have first calculated the bulk system and determined the 
lattices. The obtained values (against experimental ones 
in parentheses) are 4.09 (4.02) A for LiF, 5.13 (5.13) A 
for LiCl, 5.95 (6.00) A for Lil and 5.63 (5.63) A for NaCl. 

We have then introduced a slab model, where we put 
a jellium in the middle of the unit cell sandwiched from 
top and bottom by an alkali halide layer (Fig.l). Start- 
ing from this structure as an initial state, we have per- 
formed structure optimization. To confirm that the sim- 
plified model with the jellium capture main features in 
heterostructures with real metals, we have also studied 
the system where five layers of Cu is considered in place 
of the jellium. The lattice mismatch between LiCl and 
Cu happens to be small, so that we can take a small 
unit cell. The cut-off energy of the plane-wave expansion 
for the wave functions are 20.25 Ry for LiF, LiCl and 
Lil/jellium, 25.00 Ry for NaCl/jelhum, and 56.25 Ry for 
LiCl/Cu. The atomic configurations and the correspond- 
ing electronic states in the ground states are obtained 
with the conjugate gradient scheme^^. 



The local density of states is calculated as 
|(/)i(a;, J/, 2;)p, where (pi's are eigenfunctions and 
the summation is taken over a window Ep < Ei < 
Ep + 0.5eV with Ep being the Fermi energy. We have 
set the window above Ep to compare with the exper- 
iment, where the NEXAFS'^^ detects the unoccupied 
states. 



III. RESULTS 
A. LiCl 

Let us first discuss the case of LiCl/jellium. 
Since MIGS have been observed both in 
LiCl/Cu(commensurate) and LiCl/Ag(incommensurate)^ 
the commensurability in the lattice constants across the 
heterointerface should not be essential for the formation 
of MIGS, as we shall confirm below from the results for 
LiCl/Cu and LiCl/jellium. So we start with the jellium 
model. 

We show in Fig. 2(a) the band structure of LiCl/jellium 
where LiCl layers, each one monolayer thick, sandwiches 
the jellium that has Ts = 2.5. We have also displayed the 
band structure for an isolated LiCl monolayer. From the 
comparison we can see that the band structure of LiCl is 
not severely affected by the presence of the jellium. 

It is well-known that GGA or LDA generally underes- 
timate the band gap, while GW approximation^'^'^'* im- 
proves this. For instance, the band gap for bulk LiCl 
is 9.2 eV in GW, which is very close to the experimen- 
tal value of 9.4 eV^^. However, the band-gap underes- 
timation may be amended via the self-energy correction 



(i.e., without any corrections for LDA wave functions), 
the shape of the wave functions should be reliable even 
in LDA, hence the present calculation for LDOS around 
the Ep is expected to be a good approximation. 

In Fig.2(b,c) we have plotted the local density of states 
(LDOS) integrated along x or y-directions. We can see 
that LDOS has a significant (~ lattice constant) tail with 
a Pz character on the halogen site in addition to a usual 
evanescent- wave like part that exponentially decays away 
from the interface. In the experiment ""^^j while there is no 
information on Li sites since the element-selective NEX- 
AFS was tuned for CI in that experiment, an unoccu- 
pied state on CI sites is observed near Ep. They have 
a Pz character, and spread with a decay length of a few 
AThese properties are totally consistent with the present 
results. 



(a) LiCI/jellium (rs=2.5) 
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FIG. 2. (a) Band structure for the system where LiCl 
monolayers sandwich the jellium having Vs — 2.5 (left panel) 
and that for an isolated LiCl (i.e., rs = oo, right panel), (b) 
Local density of states (LDOS) for Ep < E < Ep + 0.5eV 
integrated over y-direction. The solid (open) circle denotes Li 
(CI), (c) LDOS integrated over both x and y directions. The 
shaded region indicates the insulator region. 



(a) LiCI/jellium(rs=2.5) 
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FIG. 3. A plot similar to Fig. 2, for a system where LiCl 
double layers sandwich the jellium having Vs = 2.5. 

If we now replace the jellium with (five layers of) Cu 
in Fig. 4, we can see that the LDOS for the LiCl/Cu has 
a long tail just as in the result for the jellium model, so 
that the characteristic feature of MIGS should not be an 
artifact of the jellium model. 



(a) LiCI/Cu 




FIG. 4. A plot similar to Fig. 2, when the jellium is replaced 
with five layers of Cu. 



To confirm that the penetration depth of the tail of 
MIGS is as large as a few A, we have also calculated a 
larger system, where LiCl double layers sandwich the jel- 
lium. From the band structure and the LDOS in Fig. 3 we 
can see that the MIGS do not penetrate into the second 
layer of LiCl. 



Why MIGS have amplitudes only on Li sites can be 
explained as follows. In general, MIGS can be divided 
into two groups: one having a conduction-band charac- 
ter and the other having a valence-band character in the 
insulator side of the interface. Noguera et aP have in- 
troduced the energy Ezcp, above which MIGS have a 
valence-band character. While the charge neutrality in 
the insulator side is satisfied when i?zcp is equal to Ep, 
charge transfer from the insulator to the metal occurs 
if Ezcp > Ep. In Fig. 5, we plot the difference in the 



valence-charge density 5p between the LiCl+jeUium and 
the isolated LiCl. We can see that, while 5p is positive 
for most of the insulator region (0 < z < 2), it is negative 
around the CI site. Thus Ezcp should be higher than Ep 
in this system, and MIGS should have a tail on the anion 
sites. 
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FIG. 5. (a) The difference in the valence charge density, 
5p = pi — p2, integrated over y-direction, where pi is the 
density for LiCl+jeUium and p2 for the isolated LiCl. (b) 5p 
integrated over both x and y directions. Inset is a blow-up 
for —0.05 < 5p < 0.05. Shaded regions indicate the insulator 
region. 



Let us move on to the dependence. In Fig. 6, we show 
the result when is increased from 2.5 to 6.0. We can 
see that, although LDOS at the interface is smaller than 
that for = 2.5, the decay length for = 6 is similar to 
that for Vs — 2.5. Thus we may expect that penetration 
depth does not sensitively depend on the nature of the 
metal. 



and Fig. 9 for NaCl. The results are similar to that for 
LiCl(Fig.2). Namely, MIGS have a long tail on anion 
sites with a penetration depth A as large as the layer 
spacing Oz = a/2. The lattice constants are a = 4.09 A 
for LiF, 5.13 A (LiCl), 5.95 A (Lil) and 5.63 A (NaCl), 
and A for Lil is indeed the largest (1.5 times as large 
as that for LiF). Note that anions with higher electron 
affinity have smaller ionic radii in general, so that alkali 
halides with largeer band gap have smaller lattice con- 
stants. 

(a) LiF/jellium(rs=2.5) 
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FIG. 7. A plot similar to Fig.2 with LiCl replaced by LiF. 



(a) LiCI/jellium(rs=6.0) 
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FIG. 6. A plot similar to Fig.2 when is increased from 
2.5 to 6.0. 



B. Other alkali halides 



(a) Lil/jellium(rs=2.5) 




FIG. 8. A plot similar to Fig.2 with LiCl replaced by Lil. 



Let us move on to the cases of other alkali halides. 
We have performed calculations for LiF, Lil and NaCl. 
The results are shown in Fig. 7 for LiF, Fig. 8 for Lil, 
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FIG. 9. A plot similar to Fig.2 with LiCl replaced by NaCl. 



IV. POSSIBILITY OF EXCITON MEDIATED 
PAIRING 

One of the most fascinating possibilities which should 
be explored for heterointerfaces is, we believe, the 
exciton-mediated superconductivity. The exciton mech- 
anism has been originally proposed by Little-'^^ for metal- 
lic, quasi-one dimensional spines of conducting electrons 
to which dielectric organic molecules are conceived to be 
attached chemically. Subsequently, two-dimensional ver- 
sion of this mechanism has been proposed^^"^*, where a 
metallic film covered by a dielectric coating was consid- 
ered. 

To realize the exciton mechanism in general, carriers 
have to coexist and strongly interact with the excitons. 
Thus the penetration length must be sufficiently large, 
which is why it has been believed that the band gap of the 
coating material should not be too large. On the other 
hand, Inkson and Anderson^® have pointed out that, if we 
adopt the model dielectric function for typical semicon- 
ductors, exciton-electron coupling becomes small, par- 
tially because the exciton has a large energy dispersion 
so that the phase space over which the mechanism works 
for the formation of paring is restricted. In fact, in Ref.^'' 
the energy dispersion of the exciton is ignored. In other 
words, they have considered a Frenkel exciton, which is 
usually unrealistic for (narrow gap) semiconductors. 

Now, the present result for MIGS at metal/insulator 
interfaces can reconcile the very dilemma. Namely, while 
a conventional wisdom tells that the penetration depth 
should be extremely small for metal/ wide- gap insulator 
interfaces, the present result suggests that the depth is 
of the order of the lattice constant for the MIGS, so that 
a substantial portion of the carriers can interact with the 
exciton in the insulator side. Thus it becomes an intrigu- 
ing problem to estimate the electron-exciton coupling for 



the alkali halide/metal system, so we do this following 
the argument of ref.^^. 

According to ref.^^, the exciton-electron coupling Acx is 
estimated as s^'^[D / L)^.[ujp/ ujg)"^ , where s is the screen- 
ing factor for the Coulomb interaction, 7 the ratio be- 
tween the LDOS at the interface and that for the bulk, 
D the average depth of penetration, L the thickness of 
the metal film, /i the screened Coulomb interaction av- 
eraged over the Fermi surface multiplied by the density 
of states at Ep, Wp the electronic plasma frequency, and 
LOg the band gap of the insulator. If we assume that the 
dielectric function can be approximated as 



(^optical ~ 1)1^0 



where ojq i^g is the typical exciton energy, and 
^optical the optical dielectric constant, then Acx becomes 
s^'^{D/L)fx{e optical - !)■ 

Let us first look into the factor ^D/L. This factor is, 
in the simplified model due to Bardeen, the fraction of 
the time over which the metal electrons spend in the in- 
sulator. In the more realistic model we have at hand, we 
can more precisely integrate the LDOS in the insulator 
region (= pi, shaded region in Fig. 2(c) etc) since MIGS 
do not possess a monotonic (exponential) LDOS, and we 
can then divide the quantity by that for the metal re- 
gion (= pm) to obtain the factor corresponding to / L. 
Namely, the factor should read b = 2pi/ pm, in place of 
^D/L, where = a/2 (see Fig.l). The present model 
is a sandwich structure where the metallic layer is cov- 
ered by the dielectric coatings from both sides, so that 
the factor 2 appears in r.h.s. Here we set that the thick- 
ness of metallic layer to be 10 A. The resulting values are 
b = 0.08 (for LiF), 0.18(LiCl), 0.16(LiI), and 0.12(NaCl). 

On the other hand, (eopticai — 1) are 0.92, 1.68, 2.4 and 
1.33 for LiF, LiCl, Lil and NaCl, respectively^^. There- 
fore, if we adopt the same values for s(= 1/2), p{= 1/3), 
7(= 1/2) as in Ref.^"^, Acx can become as large as 0.06 
for Lil, while we have Acx = O.Ol(LiF), 0.05( LiCl), and 
0.03(NaCl). 

If we adopt weak-coupling conventional superconduc- 
tors, e.g. Al, as the metal, the superconducting transition 
temperature can be roughly estimated^*'^^ as 



T, = (?cxp - 



1 



(Aph + Acx — M*) / ' 



where Aph is the electron-phonon coupling, Acx the 
electron-exciton coupling, and p* the renormalized p. So 
9 is roughly expressed as 

/duj 
(logCj) [pcx(^^) + Pph(^)] — , 

where pph (pcx) is the density of phonon (exciton) states. 
In general, pph has a peak in a low energy region, so 
that, when phonons and excitons coexist, the contribu- 
tion of the phonon part to should be dominant than 



that of the cxciton part. Thus basically scales with 
that of phonon, but the exciton contribution on top of 
that can enhance the Tc significantly. For example, when 
Aph ^ 0.4, fx* ~ 0.05, and Aex — 0.1, becomes 1.6 
times greater than that in the bulk. Thus we conclude 
that heterointerface between metal and Lil or LiCl can be 
a promising ground for detecting the existence of exciton- 
enhanced superconductivity. 

In addition, the proximity of different systems enables 
us to expect entirely different mechanism. Namely, two 
of the present authors^^ have earlier proposed that a class 
of repulsively interacting systems that consist of a carrier 
band and an insulating band can become superconduct- 
ing, where the system can effectively mapped to an at- 
tractive Hubbard model. The point is, while the conven- 
tional boson-exchange pairing arises due to a retarded at- 
traction, superconductivity in the carrier-insulator model 
occurs due to a nonretarded attraction. In their model, 
the metallic band is assumed to interact only with the 
anion sites (or to interact only with the cation sites). 
At the interface between an alkali halide and a metal, 
MIGS have their amplitudes on anion sites, so that the 
MIGS should interact primarily with anion sites, which 
just corresponds to the situation considered in Ref.^^. 
The detailed study on the combined effect of the nonre- 
tarded attraction (which decreases /i*) and the exciton 
effect (which increases Aj is an interesting future problem. 



V. SUMMARY 

We have studied the electronic properties of metal- 
induced gap states(MIGS) at interfaces between various 
alkali halides and a metal. We have performed first- 
principles density functional calculation to show that (i) 
in addition to a usual evanescent-wave component, MIGS 
have a tail on halogen sites with p^-like character, whose 
penetration depth is as large as the lattice constant of 
the alkali halide, (ii) while A does not significantly de- 
pend on the carrier density of jellium, it depends on the 
energy gap of alkali halides, with Alif < -^Lici < -^Lii- 
We have estimated the electron-exciton coupling, and 
have found that we may envisage that insulator /metal 
interfaces such as those discussed in the present study 
may provide a possible ground for exciton-mediated su- 
perconductivity. 
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